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Abstract 

Calculations for a set of nuclear multifragmentation data are made using 
a Canonical and a Grand Canonical Model. The physics assumptions are 
identical but the Canonical Model has an exact number of particles, whereas, 
the Grand Canonical Model has a varying number of particles, hence, is less 
exact. Interesting differences are found. 
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I. INTRODUCTION 



In experiments whose goals were to investigate the role of isospin in fragment yields [|TJ, 
the following interesting features have been observed. If we compare the central collisions of 
two heavy ion systems, 1 and 2, which are similar in all aspects of the reactions except for the 
neutron and proton composition, the isotope yield ratios, Y 2 (n,z)/Y!(n, z), where Yj(n,z) 
is the yield of the isotope with neutron number n, and proton number z, from reaction i, is 
found to exhibit an exponential relationship as a function of n and z |],|2| 

Y 2 (n, z)/Yi(n, z) = Cexp(a n n + a p z) (1.1) 

where C is an overall normalization constant and a n and a p are fitting parameters. This 
phenomenon is termed isoscaling, a strong evidence that the processes are statistical. 

Related to isoscaling is the exponential dependence of the mirror nuclei ratios on the bind- 
ing energy. In Figure 1, the ratios of yields of mirror nuclei: Yi(t)/Yi( 3 He), Yi( 7 Li)/Yi( 7 Be) 
and Y i ( u B)/Y i ( n C) for central collisions of 12i Sn + 124 Sn (solid points) and 112 Sn + 112 Sn 
(open points) at 50 MeV per nucleon are plotted as a function of the binding energy dif- 
ference, AEb- These ratios fall approximately on an exponential. Many statistical models 
such as the grand canonical model |3|,|j] of multifragmentation predict both the isoscaling 
and mirror-nuclei ratio dependence. 

Experimental evidence suggests that multifragmentation occurs when the heated matter 
expands to density about 1/3 of nuclear matter density H and the time scale for the emission 
of fragments is short, between 50 to 100 fm/c ||. Most successful statistical models that 
describe multifragmentation data assume a freeze-out volume at which composite yields are 
to be calculated entirely according to phase-space [0,0. If the dissociating system is very 
large, then grand canonical simplification can be employed MM- According to this model, 
the average number of composites with neutron number i and proton number j is 

< rii j >= exp(/3(i//„ + j> p ) (1.2) 

where fi n , /i p are neutron and proton chemical potentials and 

2 



Uij = ^(2vr(^ + j)mTf/ 2 (2s + l)z int exp(/3£ s ) (1.3) 

is the partition function of one composite. (3 is the inverse temperature. For mirror nuclei: 
% = k + 1, k and j = k, k + 1 we should simply have 

< n k +i,k > _ exp (^ _ ^ ) exp(/5AE jB ) (1.4) 

< ^fc,fc+i > 

Thus the log of the ratios of the yield will be linear with respect to AEb which is approxi- 
mately obeyed by data. However a more close inspection raises another issue. 

According to Eq. 1.4, one can deduce the value of /3 = 1/T from the slope of the line. 
Indeed for the lines drawn in Fig.l, the temperature T is less than 2 MeV. For such a low 
temperature, the model of simultaneous breakup model [7,8] should not be appropriate. 
In addition, such low values are in direct contradictions with temperature measurements 
obtained from isotope yield ratios. The isotope yield temperature is about 5 MeV for the 
Sn+Sn systems PJTUfl. To resolve the discrepancies between temperatures observed, it is 



necessary to explore details of the exponential behaviour of the mirror nuclei. 

II. CANONICAL VS. GRAND CANONICAL MODELS 

In recent years, the grand canonical model has been replaced by a canonical model. The 
physics assumptions are still the same but we no longer have to assume that the system is 
large. This is a technical advancement; the details have already been described in several 



places [11-13]] so we will not repeat these here. The model has been used to fit the isotope 



data [12,14]. Surprisingly, isoscaling which follows naturally from the Grand Canonical 



model, emerges also in canonical model [^4[ . In this article, we will investigate why certain 
results from the canonical model resemble those from the grand canonical model and what 
are the differences. We will also investigate the relation between the canonical temperature 
and the temperature obtained based on the simpler grand canonical rules. 

The yield of the composite which has k + 1 neutrons and k protons is given in the 
canonical model by 



QN-k-l,Z-k / n , \ 

< nk+i,k >— <^k+i,k y\ 

Wn,z 

Here N, Z refer to the number of neutrons and protons of the disintegrating system. Qn,z 
is the canonical partition function of this system. Similarly, QN~k-i,z-k is the canonical 
partition function of the residue system which has N — k — 1 neutrons and Z — k protons. 
The ratio of the yields in the canonical model is then given by 

< n k+l,k > _ ^>k+l,k x QN-k-l,Z-k ^ 2) 

< nk,k+l > ^k,k+l QN-k,Z-k-l 

The first factor leads to gxp(PAEb)- We note in passing that for mirror nuclei AEb = AEc, 
the change in Coulomb energy. If we assume a uniformly charged sphere, then AE C = 
I Roa 1/3 ^ Z ^ 2 — ^ = -^^d 2 ^ 3 MeV where a is the composite mass number. For light nuclei 
0.72a 2//3 MeV does not fit the data very well. We note for later use that 0.235a MeV fits 
AEb between a = 7 and a = 15 better. 



The exact expression for the canonical partition function Q jv,z used in |14| does not allow 
us to investigate easily the features we want to study. Since the ratios are very simple in the 
grand canonical ensemble and since there is a connection between grand canonical partition 
function Z gr (X n , X p ) and the canonical partition function Qn,z, we find it convenient to 
exploit this relation. In the present problem, the grand canonical partition function is given 
by 

Z gr (X n ,\ P )= E e^+a**" x (2.3) 

k,l,n kl i rc > 1 

The expression for logZ gr (X n , X p ) is 

logZ gr (X n , X p ) = J2 exp(fcA n + IX P ) x u k>l (2.4) 
k,l 

The canonical partition function can be obtained from Z gr by Laplace inverse: 
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Q NZ = —— / / e ^+ix n )N e -(x p +ix p )z e io g z gr (x n +ix n ^ +l x p ) d ^ nd ~ x 

(27r) 2 I - I - 
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While this expression is true for any A n and X p , the saddle-point approximation consists in 
choosing the values of X n and X p such that the kernel maximizes at X n = and X p = and 
making a Gaussian approximation for the integrand around this maximum. The result is 



Q N ,Z « e -(^N+X p Z) e logZ gr (X n ,X p ) j ^ + |^|l/2) (3 5) 

where the values of A n and A 2 are such that the average numbers of neutrons and protons 
as obtained from the grand canonical ensemble are N and Z, i.e., 

N = dlogZ gr (X n ,X p ) _ ^ ^ 

dX n 

z _ dlogZ gr (X n , X p ) ^ ^ 

dX p 

The elements of the determinant are given by: a 1} \ = 9 1 q^ 9T , Qi,2 = ^2,1 = ^°Q^ gr and 

_ d 2 logZ gr 
fl 2,2 - d 2 Xp ■ 

Eq. (2.2) now takes the form 

< T2i.ii h > -(X n (N-k-l)+X p (Z-k))+logZ gr (X n ,X p ) 

«+!,« ~ pP^Eb y (9 Q\ 

< U ktk+1 > e -(X' n (N-k)+X' p (Z-k-l))+logZ gr (X' n ,X' p ) K^J 

Here we have omitted the ratios of the determinants \det^l 2 because their effects will be 
negligible. Eq. (2.9) will reduce to the standard grand canonical result if we set X n = 
X' n ; X p = X' p and take these values from a system which has the average number of neutrons 
to be N ( rather than N — k — 1 to get X n and iV — k to get A^ ) and the average number 
of protons to be Z ( rather than Z — k to obtain X p and Z — k — 1 to obtain X' p ). 

For a better estimate, let us write A^ = A n + AA„; X' p = A p + AA P . Assuming lowest order 
expansion is valid we can get (depending upon whether we expand logZ gr (X n , X p ) in terms 
of logZ gr (X' n , X' p ) or vice versa): 

e (logZ gr (X n ,X p )-logZ gr {X' n ,X' v )) _ e -AX„(N-k-l)-AX p {Z-k) Qr e -AX„(N-k)-AX p (Z-k-l) _ 

Eq. (2.9) can be reduced to < Ttfc+1 - fc > « e^ AEB x e Xn ^ m e P AEB x e x '"~ x 'p 

H v ' <n k ,k+i> 

We will use 

< Uk+1 > k > ~ e P^E B x e (X n +X' n -X p -X' p )/2 £2.10) 

Eq. (2.10) looks just like a grand canonical result but with an important difference. In 
the usual grand canonical model A n , A p would be calculated just once, from eqs. (2.7) and 



(2.8) where iV and Z are the neutron and proton numbers of the disintegrating system. By 
contrast, A n , A p etc. of eq. (2.10) are calculated from eqs (2.7) and (2.8) for each k and the 
left hand sides of eqs. (2. 7) and (2.8) are given by iV — k — 1 and Z — k respectively. The 
quantity A n — A p etc. increases with k with the result that if we try to interpret the canonical 
results within a usual grand canonical framework one ends up with a larger (3, that is, a lower 
T. This is demonstrated in Fig. 2 where it is shown that although the temperature used for 
the canonical calculation (hence the true temperature) is 5 MeV, deducing the temperature 
from the slope of the mirror isotope yield ratios (as one would do in a grand canonical 
formalism) one would arrive at a significantly lower temperature. The best fit (solid line) to 
the calculated values from the canonical model (solid points) yield a temperature of 3.395 
MeV. In the figure we also show that the approximation of Eq.(2.10) as shown by the star 
points, works quite well. 

The dependence of A n — A p on k and N, Z where 2k + 1 is the mass number of the emitted 
particle and N, Z gives the size of the emitting system can be pinned down further. Let 
A N , Ap be the fugacities of the system N, Z. We will write A n = A N +d\^ and A p = A P +d\p. 
We then have N — k — 1 = ^2icu it j exp(iA n + j\ p ). Expressing A n , A p in terms of A^Ap 
and approximating exp(dA P ) (1 + dA P ) etc. we get — k — 1 = AdA N + BdA P where 
A and B are constants: A = J2i 2uJ ij exp(iA]y + jA P ) and B = ij^ij exp(iAjv + jA P ). 
Similarly starting from N — k = exp(iA n + j X p ) and expanding as above we get 

—k = BdA N + CdA P where C = J2 j 2 ^i,j^v{^N + jA P ). One can now express dA N) dA P 
in terms of the constants A, B and C. We get 

Q C -\- B 

\ n -\ p = A N -A P + B2 _ AC k+ B2 _ AC (2.11) 

We can do a similar analysis for A^ — \' p . Finally we get (compare eq.(2.10)) 

(A n + A; - A p - Ay/2 = A N -A P + l JL^L (2k + 1) (2.12) 

Eq.(2.12) says that the correction grows like 2k + 1 = a, the mass of the composite. The 
correction would diminish as the disintegrating system (N, Z) grows. The constants A, B 
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and C are positive definite and each will become larger and larger as the disintegrating 
system becomes larger. The correction would disappear in the thermodynamic limit. The 
actual values of the constants A, B and C for a finite system depend on many factors: the 
symmetry energy, the Coulomb energies and N, Z of the disintegrating system. 



III. THE ALBERGO TEMPERATURE 

The Albergo formula 0] has often been used to extract a temperature from experimental 
data. The formula is exact if the following two assumptions are valid: (1) the populations of 
various states are given by the grand canonical model and (2) the secondary decays which 
will alter these primary populations can be neglected. Define a ratio R 

y(A,zQ/y(A + i,z») , * 

Y{A ji Z J )/Y{A i + \ i Z j ) 1 • 1 

where the Y's are the yields in the ground state. Then, the temperature is given by 

T -^W) < 3 - 2 > 

where B is related to binding energies and s to the ground state spins: 

B = BE(Ai, Zi) - BE(Ai + l,Z t )- BE(Aj, ZA + BE(Aj + 1, Z 3 ) (3.3) 

[2g(A 3 -,Z 3 -) + l]/[2g(A 3 - + l > Z 3 -) + l] 

[2S(A, Zi) + l]/[2S(Ai + 1, Zt) + 1] 1 ' 1 

Even if the grand canonical model is exact, the change of populations due to secondary 
decays can cause eq.(3.2) to give significantly different temperatures from the true grand 



canonical temperature. This was studied in detail in [15[]. It was shown that for large values 
of B (eq.(3.3)), the difference between apparent temperature and the true grand canonical 
temperature decreases. This suggests that while using the Albergo formula to deduce a 
temperature from experimental data, it is advisable to use pairs that will lead to a large 
value of B. 
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Our objective here is different and is complimentary to the study made in [15|. The 
canonical model is obviously more rigorous than the grand canonical model. However, if 
canonical values for R are used, eq. (3.2) is no longer strictly correct. Using the primary 
yields, we explore the differences between the deduced temperatures from eq.(3.2) compared 
to the actual temperature used in a canonical model. This is shown in Fig. 3, we find 
that the errors decrease with increasing B. The inset in Figure 3 shows the deviation 
of the canonical Albergo temperature for B greater than 10 MeV. Most of the predicted 
temperatures are slightly lower than the actual temperature of 5 MeV. The deviations arise 
from the differences between isotope yields predicted by the canonical and grand canonical 
models Not surprisingly, the conclusions of |15j can be applied here. 



IV. THE SCALING LAW 



The last quantity we want to investigate is a ratio of two ratios: < n '+ m - fc > 2 4- <ni - k > 2 and 
see if this falls on an exponential as in the grand canonical ensemble. Here the subscripts 1 
and 2 refer to two systems : ( for example: 2 refers to central collisions of 124 Sn+ 124 Sn and 1 
to central collisions of 112 Sn+ 112 Sn at 50 MeV/A energy). As this involves two ratios and two 
different systems, the analysis is considerably more complicated than what we considered 
before. The ratio R we are after is given by 

,, (2) n n 

D u l,k VJV2-!-m,Z2-i <°iNl-l,Zl-k /a i\ 

K= — (2) (T) n T\ t 4 --U 

u lk ^l+mk i *N2-l,Z2-k WNl-l-m t Zl-k 

For central collisions at the same beam energy per particle, the u factors will give unity. 
Employing the saddle-point approximation and setting the ratios of the det's as unity as 
before, we can consider 

Q ~(X n (N2-l~m)+X p (Z2~k))+logZ gr (\ n ,X p ) 
N2— I— m,Z2— k c 



QN2-1 Z2-k ' e -{\' n {N2-l)+\' p {Z2-k))+logZ gr {\' n ,\' p ) 



(4.2) 



Similarly, 



n -(\ n {Nl-l-m)+\ p {Zl-k))+logZ gr {\ n ~\ p ) 

Qm-i,zi-k ~ e -(K(Ni-i)+~y p (zi-k))+iogz gr Cx n ry p ) ^ ' ' 



We can now indicate how the grand canonical results are recovered. We set \ n = \' n ;\ n = X n 
and A n — A n = AA„ then the ratio achieves the exponential character: R = exp(mAA„). 
Experimentally it is found that the relationship R = exp(am) where a is a constant inde- 
pendent of I and k is quite well respected. This is not so obvious from eqs. (4.1), (4.2) and 
(4.3). We are therefore required to investigate this near independence of the constant a. 

If we write in Eq.(4.2) A^ = A n +AA n and expand logZ gr (X n , X ) in terms of logZ gr (\ n , X p ) 
and keep lowest order corrections, the right hand side of eq.(4.2) is simply e mX ' n . In a 
similar fashion, the right hand side of eq.(4.3) is e mX ' n , so that the ratio R of eq.(4.1) is 
exp(m(A^ — X n )) where, of course, the values of A^, A^ are chosen to give neutron numbers 
N2 — I and Nl — I and proton numbers Z2 — k and Zl — k respectively. Our next task is 
to verify that X n — X n is approximately independent of I and k. 

We have four equations: 

]T ie iX '" +jX 'vuj iyj = N2-l (4.4) 
je^+'^Wij = Z2 - k (4.5) 
J2ie fx '" +j ^uj tyj = Nl -I (4.6) 

Y,je*> +j %u i>j = Zl-k (4.7) 

Let X n = X n + 5X n and A^ = A^ + 5X p . From Eqs. (4.5) and (4.7), retaining terms to lowest 
order in SX P and 5X n we obtain 

$K E fe^+^Uij + 5X n E ijJ^+^Uij = Z2-Z\ (4.8) 

In a similar fashion from Eqs. (4.4) and (4.6) we can obtain 

5X n i 2 e lX '" +jX 'vuJ hJ +5X P Y,i 3e fx '" +f ^uj hj = N2 - Nl (4.9) 

Eqs (4.8) and (4.9) can be solved for SX n and 5X P and in the lowest order these value are 
independent of / and k but depend upon N2,Z2, Nl and Zl. To this order R of eq.(4.1) 
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is independent of I and k as it is in the usual grand canonical ensemble. This is seen to be 
obeyed in experiments to a good approximation. 

Instead of eqs.(4.4) to (4.9), one may also consider the following approximation trying 
to relate to the grand canonical ensemble. Recall that A^, A^ are the fugacities of a system 
which has N2 — I neutrons and Z2 — k protons. If we denote the fugacities of the system 
which has N2 neutrons and Z2 protons by Aat 2 ,A^ 2 and employ the same approximate 
methods used in the discussion leading to eq.(2.11), we get A^ = Aat 2 + ^-mc-i Similarly 
X' n , \' p refers to a system which has Nl — I neutrons and Zl — k protons. In an obvious 
notation we also get A^ = Ajvi + Jrf^jf^ ■ The quantity of interest is 



Because of cancellations in the above equation, results again approximate the grand canon- 
ical result quite closely. 



V. SUMMARY 

In summary, we have explored several experimental observables which are sensitive to 
the isospin effects in multifragmentation. We find that the mirror ratios, isoscaling and 
temperatures calculated in canonical model behave similarly as those predicted with the 
grand canonical model with one significant difference: the temperature deduced from the 
calculated observables with the canonical model using the rules based on the grand canonical 
model can be significantly different from the true temperatures. 
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FIGURES 

FIG. 1. Isobar ratios for three pairs of mirror nuclei obtained from the central collisions of 
124 Sn+ 124 Sn (solid points) and 112 Sn+ 112 Sn (open points) collisions. The lines are best fit of Eq. 
1.4. 

FIG. 2. Exact Canonical Model calculations for a system of neutron number N=104 and Z=70 
using a freeze-out density of one-quarter of normal density. This simulates central 124 Sn+ 124 Sn col- 
lisions. Lower values of N and Z used in this calculation reflect effects of pre-equilibrium emissions. 
The ratios of yields of mirror nuclei are plotted for a = 1, 3, 7, 9, 11 and 13. The results for 3,7 and 
11 can be compared with the experimental results (Fig.l). T actua i=5 MeV is the temperature used 
in the canonical model calculation; Tbestfit would be the temperature deduced if one fit the solid 
points obtained from the canonical calculations, using the grand canonical formula, Eq. (1.4). The 
results from a saddle-point Eq. (2.10) approximation are also shown. 

FIG. 3. The inverse Albergo temperature, Eq. (3.2), from the canonical model is plotted as 
a function of the binding energy difference, B. The inset shows the predicted temperature in an 
expanded scale. The dash line at T=5 MeV is the input temperature to the calculation. 
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